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We study the production of dileptons in relativistic nucleon-nucleon collisions. Additionally to the 
traditional dilepton production channels (vector meson decays, meson and A(1232) Dalitz decays) 
we included in our model as new dilepton sources the Dalitz decay of higher unflavored baryon 
resonances with spin<5/2 and mass<2.25 GeV/c^. The contributions of these new channels are 
estimated using experimental information about the decays of the resonances and have large 
uncertainties. The obtained dilepton spectra are compared to the experimental data by the DLS 
collaboration. Predictions for the HADES detector (SIS, GSI) are also discussed. In spite of the large 
uncertainties of the higher resonance Dalitz decay contributions we are able to draw the conclusion 
that these contributions are negligible compared to the other dilepton sources and do not influence 
the detectability of the <j> and lo vector meson peaks. 

PACS numbers: 13.75.Cs; 13.30.Ce; 13.40.Hq; 14.20.Gk 



I. INTRODUCTION 

The main objective of the physics program at SIS (GSI, Darmstadt) is the study of hot and dense hadronic matter 
formed in nucleus-nucleus collisions. The clearest probes for this investigation are dileptons, which leave the interaction 
volume unaffected by strong final state interactions. The purpose of the HADES detector is to study the dilepton 
spectrum produced in nucleus- nucleus, nucleon-nucleus, nucleon-nucleon, and pion-nucleus collisions. 

Theoretical investigations predict a partial restoration of chiral symmetry in hot and dense nuclear matter, which 
results in a modification of particle properties. For an experimental study of these phenomena the light vector mesons 
p, Lu, and (f) seem to be a suitable probe. The p - owing to its short lifetime - decays inside the hot and dense hadronic 
matter. Predicted in-medium modifications of lu and (f> may result in a considerable growth of their decay widths and 
there is hope that they also decay inside the hot and dense phase. On the other hand, p, uj, and (j) decay directly to 
dileptons, which allows a clear experimental study. 

Dilcptonic decays of vector mesons contribute to the dilepton invariant mass spectrum as peaks or - in the case of 
the /5 - as a wider distribution. In-medium modification of vector meson masses and decay widths may be identified 
in the dilepton spectrum as a shift and widening of the corresponding peaks, respectively. 

To identify the vector mesons in the dilepton spectrum a careful analysis of the backgrounds is needed. One of 
the possible sources of background is the Dalitz-decay of baryon resonances, TV* — > Ne~^e^ . These processes may be 
interesting also for themselves since baryon resonance properties may also be subject to in-medium modifications. In 
earlier studies only the A(1232) and iV(1440) Dalitz-decays have been considered in model calculations (the latter 
giving a negligible contribution). On the other hand, some of the higher resonances have a relatively high photonic 
decay width and are, therefore, expected to have a non-negligible contribution to the dilepton spectrum through their 
Dalitz-decay. The kinematical upper limit for the mass of dileptons resulting from the Dalitz-decay of resonances 
rises with increasing resonance mass, therefore higher mass resonances may contribute to the dilepton spectrum in a 
higher dilepton mass range where the contribution of the A(1232) is already negligible. It is particularly interesting 
to see if some of the higher mass resonances contribute essentially to the dilepton invariant mass spectrum in the (j) 
mass region since this may influence the possible identification of the (f) meson in the dilepton spectrum. Some of the 
higher mass resonances have also higher spins. In studying their electromagnetic transitions special care should be 
taken of their higher spin nature. 

Dilepton production in proton-proton collisions with beam energies between 1 and 5 GeV has been studied experi- 
mentally by the DLS collaboration , supplying an opportunity to test models of dilepton production in elementary 
processes. These models may in turn be used as an input in transport models of heavy ion collisions. 

In this paper we discuss the Dalitz-decay of baryon resonances up to spin-5/2. We inserted our results into a Monte 
Carlo model of relativistic proton-proton collisions and obtained the dilepton invariant mass spectra. 
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The lack of experimental information on the Dalitz-decay of baryon resonances forces one to use information on the 
- often also poorly known - radiative decays, N* — > iV7. This brings in several assumptions and leads to uncertainties 
of the results. The resulting uncertainty of the total dilepton invariant mass spectrum is moderate and - as we will 
see - allows one to draw the conclusion that the detectability of the <j) and uj meson peaks is not influenced by the 
baryon resonance Dalitz-decays. 

In an earlier paper we discussed the calculation of the Dalitz-decay of baryon resonances and made a comparison 
of the contributions of the various possible decay amplitudes. In Sec. ^ of the present paper we briefly review those 
results. In Sec. IIII ^ we give a short description of our Monte Carlo model for proton-proton collisions, in Sec. IIII 51 
we discuss the dilepton production channels included in the model. In Sees. IIII Ul and llil Dl we compare our dilepton 
spectra with the DLS data and give predictions for the HADES detector. Throughout the paper we set c = Ti — 1. 



II. DALITZ-DECAY OF A SPIN-J RESONANCE 



A. Kinematics and matrix elements 



Let p* , m* , and A* denote the four- momentum, mass, and helicity of the decaying resonance, respectively, and 
p, m, and A the corresponding quantities for the nucleon. Let us introduce the notation q = p* — p for the photon 
four-momentum, — for the virtual photon mass {— dilepton invariant mass) squared, and P = [p* +p)/2. 

The differential width of the Dalitz-decay of a particle is related to its photonic decay width to a virtual photon 
Tn'-.n^{M) by i 

dM^ 3^M^^^*^^^^^^)- 

Here T M-^N-^iM) can be expressed in terms of the photonic decay matrix element (A^7|r|iV*) as 



pol 

where 

(iV7|T|A^*) = -e^(iV|J^|A^*), (3) 

with the photon polarization vector, and the electromagnetic current operator. In ^ ripoi^N* is the number of 
polarization states of the N* resonance and A(a, 5, c) = a'^ + b'^ + — 2{ab + bc + ac) is the usual kinematical factor. 

Next we determine the matrix element (A^| J^JiV*) for an N* resonance of spin- J. A spin-J (J > 3/2) fermion can be 
described by a spinor-tensor field vfi" (n = J — 1/2), therefore in momentum-space matrix elements the particle is 
represented by a spinor-tensor amplitude u^^ ' ' {p* , A* ) , which has to fulfill the generalized Rarita-Schwinger relations 
in order to eliminate the different spin degrees of freedom. The electromagnetic current matrix element can be written 
generally as 

{N\J^\N*)=u{p,X)r^p,...p^uP^-P"{p*,X*), (4) 

where the form of rppj...p,^ is restricted by the conservation of electric charge, g'^(iV| A^*) =0. As a result of the 
Dirac equations u{p, \){^ — m) — and {^* — m*)^''^ '''" (p*. A*) = and the generahzed Rarita-Schwinger relations 
T/jpi - pn contains only three independent terms: 
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with 



Xpp 


= 7p9p- 


^3pp, 




Xpp 


P^,qp- 


(P ■ q)9t^p, 




App 


= q^iqp - 


M^g^^p, 


(6) 



and G = 1 or 75 for resonances with positive or negative normalities, respectively. The normality of a spin-J baryon 
is by definition P(— 1)'^"^/^ with P the intrinsic parity. fi{M^), i = 1, 2, 3 are three independent form factors, which 
should be determined from experimental data. 
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B. Parameter fitting 

In our model we neglect the dependence of the form factors, therefore fi become constant couplings. (Vector 
meson dominance does not apply here since we include the p dileptonic decay separately.) We redefine the xlip in our 
matrix elements dividing them by the appropriate power of the nucleon mass m in order to render the form factors 
fi dimensionless. We, further, write fi = egi, with e the elementary charge related to the fine structure constant 
by — Ana. We will use the real photonic decay width of the resonances to determine the dimensionless coupling 
constants gi. 

A detailed investigation of the Dalitz-decay width calculated from the 3 independent matrix elements and the 
interference terms leads to the following conclusions (see for details): 

• Those terms of the Dalitz-decay width that contain x^^p give zero contribution at M = (real photons) , therefore 
information on the real photonic width of the resonances can not be used to fix the value of the g^, coupling. 

• The remaining terms (oc gf, g| and (7152) give very similar contributions to the Dalitz-decay width. Thus an 
arbitrary choice between the xltp E^nd xftp terms causes only a minor uncertainty in the dilepton invariant mass 
spectrum, that is definitely smaller than the uncertainty arising from the poor knowledge of the photonic width 
in the case of most of the resonances. 

• The resonance mass dependence of the Dalitz-decay width may influence the resulting dilepton spectrum because 
the masses of the resonances are generated according to a Brcit-Wigner distribution in our model (see Sec. 1111 A|l . 
but we use a single coupling fitted to the photonic width of a resonance with its mass equal to the peak value. 
However, the Breit-Wigner distribution ensures that resonance masses stay close to the peak value and hence 
the uncertainty of the dilepton spectrum caused by the different resonance mass dependence of the xjip s-^d xfip 
is small. 

These results encourage us to make the assumption that g2 — gs — 0, i.e., only Xpp (the term containing the lowest 
power of external momenta) contributes to the Dalitz-decay of baryon resonances. 

To demonstrate the relevance of the spin-parity of baryon resonances in their Dalitz-decays we show in Fig. ^ the 
differential Dalitz-decay width of hypothetical resonances with the same mass and photonic width but with different 
spin-parities. It can be seen from the plot that at the high dilepton mass end of the spectrum spin-parity 1/2+ and 
3/2— resonances give the largest while 5/2— give the smallest contribution to the dilepton spectrum. 

III. DILEPTON PRODUCTION IN PROTON-PROTON COLLISIONS 

A. Tiie model 

We use a resonance model to describe particle production in proton- proton collisions: baryon resonances are pro- 
duced in the first step, which then decay - possibly in a multi-step process - into the final state particles. The 
resonances are treated as on-shell particles in their production processes, but their mass is generated according to a 
Breit-Wigner distribution. We included in our model all unflavored baryon resonances with mass<2.25 GeV, and with 
a status of at least 3 stars according to the Review of Particle Physics i.e., A(1232), iV(1440), iV(1520), iV(1535), 
A(1600), A(1620), 7V(1650), iV(1675), 7V(1680), 7V(1700), A(1700), iV(1710), iV(1720), A(1905), A(1910), A(1920), 
A(1930), A(1950), Ar(2190), iV(2220), 7V(2250), and some additional resonances with weaker status, namely A(1900), 
Af(2000), and 7V(2080). In our model the following decay channels of baryon resonances are taken into account: Ntt, 
Nt], Nrnr (including Np, A(1232)7r, 7V(1440)7r, and 7V"cr"), Nuj, and KK. 

We fitted the resonance properties (branching ratios, full widths, and masses within the errors given in the Review of 
Particle Physics |^] ) to experimental data of pion induced reactions. Cross sections of resonance production have been 
fitted to one-pion, two-pion, 77, uj, p, and kaon production in proton-proton and proton-neutron collisions. Formulae 
for resonance and meson production cross sections used in the fits are essentially the same as in the BUU model 
established by Gy. Wolf (see 5] and references therein). We needed to include in the model the resonances with 
weaker status mentioned above in order to obtain a good fit. For further details of parameter fitting see 

B. Dilepton production 

We used the above model to calculate the dilepton invariant mass spectrum of proton-proton collisions. The dilepton 
sources considered are the direct dilepton decay of p''\ oj, and 4> vector mesons and Dalitz-decay of 7r£, rj, and w mesons, 
and baryon resonances. For the meson Dalitz-decay channels we used the same expressions as in [5j. 
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The differential cross section of dilepton production in the direct decay of a vector meson is given by 



dM dM 



BRy_^,^,-{M). (7) 



Here BRy^f.+ e- {M) is the dilcptonic branching ratio of the vector meson V , which is approximated in the case of 
the p as 

i?i?p^e+e- {M) = r^^^^(Af) + r|!J+l^(M) + r,^.^(M)' 

where we adopt the mass dependence of the pionic and dileptonic p width from 7] while we keep T p^-j^^ constant, 
determining its value from the Breit-Wigner with and branching ratio listed in the Review of Particle Physics 0. 
The appearance of Fp^g+g- and Fp-,^^ in the denominator of JSJ is only relevant close to the two pion threshold 
M = 2m^, where Fp_,^7r becomes small. 

We keep the dileptonic branching ratio of the narrow oj and (j) resonances constant and use their PDG values. </> meson 
production is not handled by our resonance model and the production cross section is not well known experimentally. 
This leads us to adopt here the (f) meson production cross section calculated from a one-boson exchange model Q . 

We calculate the dilepton spectra from the Dalitz-decay of those of the above mentioned baryon resonances that 
have spin<5/2 and are known to couple to the iV7 channel according to the PDG, i.e., A(1232), A^(1440), A^(1520), 
iV(1535), A(1600), A(1620), 7V(1650), A^(1675), 7V(1680), A^(1700), A(1700), iV(1710), iV(1720), A(1905), A(1910), 
A(1930). For resonances with spin>3/2 we use the results of the previous section. The Dalitz-decay of the iV(1440) 
resonance has already been included in Ref. Q. We treated all the other spin-1/2 resonances in the same fashion. 

To fix the electromagnetic coupling constants gi defined in Sec. ^] the photonic width of the resonances is needed. 
This we calculate using the full width of the resonances (as obtained in our resonance model) and the photonic 
branching ratios, that we take from the Review of Particle Physics Q. We made calculations using both the minimal 
and the maximal values of the photonic branching ratios given in Q . This provides us with an estimate of the resulting 
uncertainty of the dilepton production cross sections. 

Carrying out the parameter fitting described in Sec. IIII Al we arrived at two different parameter sets describing 
the experimental data with nearly the same accuracy. We calculated the dilepton spectrum with both parameter 
sets. From the point of view of dilepton production the main difference between the two parameter sets is in the 
baryon resonance Dalitz-decay channels: using the parameters of our best fit two of the higher resonance channels, 
the A^(1520) and A'^(1535) give considerable contributions to the dilepton spectrum. The contribution of these two 
channels exceeds that of the A(1232) at the large dilepton mass end of the spectrum for large enough beam energies. 
In the case of the other parameter set the A^(1520) gives a smaller contribution and its role is partially taken over by 
the A(1700), which gives a larger contribution for large dilepton masses. We used the parameters of our best fit in 
obtaining all numerical results presented in the rest of the paper. 



C. Comparison with the DLS data 

In Fig.Elwe plot the dilepton invariant mass spectra da /dM for pp collisions calculated from our model in comparison 
to the DLS data p]. We have taken into account the DLS experimental filter and the finite dilepton mass resolution 
of the detector. In the plots contributions of the following sources are shown: direct dileptonic decay of p*^, lo, and 
(f) vector mesons, Dalitz-decay of 7r°, 77, and oj mesons, and Dalitz-decay of A(1232), iV(1520) and A^(1535) baryon 
resonances. The total cross section is a sum of all sources including the Dalitz-decay of the other baryon resonances. 

A first glance at Fig. [3 shows that although the resonance Dalitz-decay channels, esp. the A^(1535) channel has a 
large uncertainty arising from the poor knowledge of the photonic branching ratios, the resulting relative uncertainty 
of the full dilepton production cross section is not dramatic. 

The dominant dilepton sources at the two lowest beam energies are the 7r° and A(1232) Dalitz-decay, and direct 
dilepton decay of p^ . At 1.04 GeV the experimental data are underestimated by the model. Comparison with results 
of other models successfully describing the data 0, ITol ITlj suggests that the reason is probably the underestimation 
of the A(1232) channel at this energy. 

For beam energies i?b > 1.61 GeV the rj channel gives the dominant contribution around 0.2 GeV dilepton mass. 
Figure |31 shows a comparison of the r] meson production cross section deduced from our model to the available 
experimental data 0,^3- Our model underestimates the two data points at Eb = 2.02 GeV {piab = 2.81 GeV). As 
a result we also underestimate the dilepton cross section in the M — 0.2 - 0.4 GeV range for the beam energies i?f, = 
1.85 GeV and 2.09 GeV (and, to a lesser extent, for Eb = 1.61 GeV). 
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The peak appearing in the dilepton spectrum of the p channel at M 0.4 GeV is a result of two reasons. First, 
low mass p mesons are produced with a high cross section via the decay of the A^(1520) resonance, therefore also 
the p mass spectrum has a peak in the same region. Second, close to the two pion threshold the width of the decay 
p —>■ TTTT becomes small, consequently the dileptonic branching ratio BRp^f,+ ^-- (M), that enters the expression of the 
dileptonic differential cross section lO, rises rapidly as M approaches 2m^. 

In the M > 0.4 GeV region the experimental dilepton spectra are more or less well described by the model for 
the beam energies 1.27 GeV < Et < 2.09 GeV with the exception of a discrepancy around M = 0.6 GeV for Eb = 
2.09 GeV. The same underestimation of the data was seen using the model ^3 and a strong pp bremsstrahlung 
contribution was named as a possible explanation. 

At 4.88 GeV beam energy channels with additional pions in the final state are known to dominate the dilepton 
spectrum [tL llollTll |. These models use a string fragmentation model to describe the DLS data. Our model, since it 
does not handle final states with additional pions, looses its validity at this beam energy. Consequently, our results 
are an order of magnitude below the experimental data. 

D. Predictions for HADES 

Figure 0] shows the dilepton invariant mass spectra calculated from our model for the beam energies Ef, — i GeV 
and 3.5 GeV. Here no experimental filter has been used. In choosing the beam energies we have taken into account 
that here we are primarily interested in the Dalitz-decay channels of higher baryon resonances, which may have an 
important contribution in the higher dilepton mass region where the contribution of the A(1232) is already small, 
i.e., around M = m^. Therefore, we have chosen beam energies above the (f) meson production threshold of Eb = 2.59 
GeV. On the other hand, one must not go beyond the beam energy where the string channel opens and the model 
looses its validity [lol |. 

In Fig. 2] we see that the dominant dilepton sources at these energies are the 7r° and t] Dalitz-decay and the direct 
dilepton decays of vector mesons. Both the N{1535) and the iV(1520) channels exceed the A(1232) contribution in 
the (j) mass region, but they remain well below the (p meson contribution. Consequently, a pronounced peak appears 
at the (f> mass in the total dileptonic differential cross section. A similar peak is seen at the w mass. These peaks 
should be clearly identified by a high mass resolution dilepton spectrometer such as HADES. We note that we found 
the same qualitative features in the case of the second resonance parameter set, namely, the (f) and uj peaks are clearly 
visible in spite of the fact that the A(1700) channel exceeds the A(1232) contribution in the (p mass region. 

IV. SUMMARY 

In this paper we investigated the production of dileptons in rclativistic proton-proton collisions using a Monte 
Carlo model. We introduced as new channels the Dalitz-decay of a number of higher baryon resonances. The lack 
of experimental information about the Dalitz decay of higher resonances causes a large uncertainty of the calculated 
dilepton spectra of these channels. We have given an estimate of the uncertainties arising from the poor knowledge of 
the iV7 branching ratios of the baryon resonances. Despite the large uncertainties we may conclude that the Dalitz 
decay of higher baryon resonances gives a negligible contribution to the dilepton spectrum. The dominant dilepton 
sources are the 7r° and 77 Dalitz-decay and the direct dilepton decays of vector mesons. 

The obtained dilepton spectra compare well with the experimental data of the DLS collaboration for beam energies 
Eb < 2.09 GeV with some exceptions. The most serious disagreement between our results and the experimental data 
can be explained by comparing the dilepton and the 77 production data at the same energies. We have also given 
predictions for the HADES detector for beam energies above the ((> meson production threshold. Our model predicts 
clear peaks at the uj and (j) meson masses, which should be easily identified by the HADES detector. 

Our model can also be used as an input to a transport model of heavy ion collisions, where in-medium modifications 
of vector mesons and baryon resonances can be studied. 
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FIG. 1: Differential width of the Dalitz-decay of hypothetical baryon resonances with the same mass (1.5 GeV) and photonic 
width (0.6 MeV), but with different spin-parities. Curves for 1/2— and 3/2+ resonances nearly overlap, similarly 1/2+ and 
3/2— curves. 



8 




FIG. 2: The dilcpton invariant mass spectra for pp collisions calculated from our model taking into account the DLS filter and 
the finite mass resolution in comparison to the DLS data. For the meaning of the various line styles see the explanation in the 
plot for f .04 GeV beam energy. The solid line represents the sum of all chaimels including those that are not indicated in the 
plots (i.e. Dalitz-decay of the other baryon resonances). In the case of baryon resonance Dalitz-decays the lower (upper) edge 
of the shaded region corresponds to the dilepton production cross section calculated using the minimal (maximal) value of the 
photonic branching ratio given in the Review of Particle Physics. The resulting uncertainty of the total dilepton spectrum is 
indicated similarly. 
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FIG. 4: Dilepton invariant mass spectrum from proton-proton collisions at various beam energies. No experimental filter has 
been used. The uncertainties of the baryon resonance Dalitz-decay channels and the total dilepton spectrum are indicated 
similarly to Fig.|21 



